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ABSTRACT 

The viscoelastic properties of dilute aqueous solutions of 
polyethylene oxide were studied using a modified simple elastometer. 
The dynamic viscosity and dynamic rigidity are calculated from the 
measurements in the frequency range from 1 to 11 Hertz for concentra- 
tions of .25%, .5%, 1%, and 1.5% of a polymer having a molecular 
weight of about 4 million. 

The results indicate that the solutions may be represented as 


a simple Maxwellian element in the frequency range considered. 
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1. Sineroduct ion 
Polyethylene oxide, a high molecular weight, long-chain linear 


molecule with the structure oe below, \ 
x Oo / AN G 
ies om coo sett 
~ 
¢ ae Om 
S ff 


has been the subject of many recent investigations. When dissolved in 


> 


water, it forms a viscoelastic solution which exhibits some of the 
characteristics of both a solid and a liquid. That is, it has an 
elastic modulus and a coefficient of viscosity. Polyethylene oxide 
solutions also exhibit a decrease in the apparent viscosity with in- 
ereasing rates of shear, which is a characteristic of so called pseudo-~ 
plastic or non-Newtonian fluids. Philippoff first found this phenomenon 
for certain fluids; he also found that the viscosity of polymers under 
oscillating deformations depends only on the frequency and not on the 
amplitude of deformation. (7) 

Non-Newtonian fluids are best characterized by a complex vis- 
cosity or impedance function y= v- Qa, Where 7a is the usual 

4 

dynamic viscosity and ¥e3= = where a is the dynamic rigidity and 
G&G) is the angular frequency of oscillation. 

The real and imaginary parts of veld give good insight into the 
properties of the solution under periodic deformation. If Ye is 

é 

zero, the fluid is perfectly viscous, and if Yes. is zero, the fluid 
is perfectly elastic. The use of a complex impedance, however, is 
useful only for oscillatory motion. 

Polyethylene oxide solutions have been the subject of two thesis 
studies at this institution. Lt. Chester found, in 1964, that these 


solutions behave essentially like water for sound waves at relatively 


high frequencies (10 to 26 megaHertz). (1) The long-chain molecules in 
solution contribute no measureable effects to the absorption or speed 
of sound. He did find, however, that shear waves in the megaHertz 
range are strongly affected by the addition of the polymer. In 1965, 
Lt. Kinnier and Lt. Reister found, using a torsional oscillator, that 
from 118 Hertz to 675 Hertz, viscosity changes with frequency, and that 
the long-chain molecules have a great effect on the properties of the 
solution. (4) 

The object of this investigation was to determine the viscoelastic 
properties of the solutions in the low oscillatory frequency range from 
1 to 10 Hertz. The polyethylene oxide used in this investigation was 
purchased from the Union Carbide Company. It is named Polyox WSR-301; 
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its average molecular weight is approximately 4 X 10° and the distri- 


bution is unknown. 


2. Equipment 

In order to measure the complex viscosity at low frequency in 
Polyox, a modified simple elastometer was constructed, following the 
designs of Goldberg and Sandvik, (3) and Markovitz, Yavorsky, Harper, 
Zapas and DeWitt. (6) Goldberg and Sandvik designed and constructed a 
simple elastometer which provided a quick and simple method for deter- 
mining viscoelastic properties of gels and soap solutions. Its opera- 
tion depends on the motion of a cylindrical bob suspended in fluid 
undergoing periodic angular oscillation. Its operation, however, is 
limited to the frequency of natural resonance for the suspended cyl- 
inder. Markovitz, et al took the basic design of Goldberg and Sandvik 
and extended its Frequency range. 

The equipment used in this investigation is similar to the elasto- 
meter of Markovitz, but utilizes an optical method for obtaining the 
needed measurements. It consists essentially of a circular cylinder 
containing a liquid, in this case Polyox, in which another smaller 
cylinder is suspended concentrically. The outer cylinder is made to 
oscillate mechanically through a small angle with a frequency that can 
be varied. The motion of the inner cylinder with respect to the outer 
cylinder can be used to determine the complex viscosity of the fluid. 

The apparatus (Figures 1 and 2) consists of a steel table with 
three leveling legs, on which a turntable with a lathe-turned shaft is 
mounted. The turntable is supported by ball bearings and is rigidly 
connected to a light aluminum rocker arm. A smooth rubber-covered 
wheel mounted at the end of the rocker arm is held by a spring against 
an eccentric cam. The eccentric cam is driven by a shaft and pulley 


system coupled to a variable speed electric motor. Since the drive 


motor is more stable in operation at high speeds, two sets of pulleys 
are used to reduce the shaft speed. The pulley system also isolates 
the steel table from any motor vibration. 

If the eccentricity of the cam is small compared to the length of 
the rocker arm, then the motion of the turntable can be considered 
sinusoidal. 

A metal can is mounted on the turntable with a smooth glass cyl- 
inder cut from a length of large diameter glass tubing fixed in the 
center. This glass tube serves as the outer cylinder and the metal can 
aS a container for a constant temperature bath. 

Several different diameter inner cylinders were constructed using 
machined and polished Lucite rod. The inner cylinder is suspended by 
a thin (No. 38) copper wire fixed to a screwjack device which permits 
the suspension height to be varied. The screwjack device can also be 
moved laterally to permit centering of the inner cylinder. 

Optical means are employed to measure the motion of the inner and 
outer cylinders. This system functions as follows: a mirror is mounted 
on the side of the outer cylinder and another mirror is mounted on the 
inner cylinder. These two mirrors are illuminated with a highly colli- 
mated beam of light from a common light source. The reflected beams 
of light are displayed on a curved screen placed three meters from the 
inner cylinder. The curved screen eliminates any distortion of the 
end points of the arcs of reflected light. Although the reflecting 
mirrors are not the same distance from the screen, the distance of sep- 
aration is small compared to the distance to the screen so that the 
error caused by this is negligible. 


Two photocells are mounted behind slits on the curved screen. 
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They can be positioned in the midpoint of the arc of reflected light 
from the oscillating cylinders. The output of each one of these cells 
is amplified and displayed on separate traces of a dual beam cathode 
ray oscilloscope (Figure 3). This permits measurement of the phase 
angle between the displacements of the inner and outer cylinder by 
using a vernier delay dial to align the output traces. The amplitude 
of the oscillation is measured directly on the face of the screen by 
measuring the length of the light arc. 

There was a slight amplitude modulation of the magnitude of both 
reflected arcs. This may have been caused by a minor wobble in the 
turntable, resulting from the light load on the supporting thrust 
bearings. 

This apparatus is useful in the frequency range from 1 to 11 
Hertz. At higher frequencies the rocker arm tends to float or vibrate 
on the eccentric cam, and at lower frequencies, the size of the inner 
cylinder, the restoring force of the wire, and their relation to the 


viscosity of the fluid becomes very critical. 
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3. Theory of Measurement 

The measuring device (Fig. 4) consists of a lucite cylinder of 
radius r suspended to depth h in a liquid contained in another cylinder 
of radius R. The outer cylinder is caused to execute simple harmonic 
oscillations through a small fixed angular amplitude. The relationship 
between the motion of the inner cylinder and the viscoelastic proper-~ 
ties of the liquid was developed by Markovitz. (5,6) 

Consider a cylinder lamina of the liquid of height h with inner 
radius z and outer radius z +dz. The torque which acts on the cylin- 
drical surface of the cylindrical lamina at z is Pe Lhaee, 
where @ is the angular velocity of the liquid at gz and Ae is the vis- 


cosity coefficient, which may be either real or complex. At z + dz, 


the torque on the surface is ene v ns) +f0 fema'fen(as) } } dz , and the 
at 


net torque on the lamina is 


ee is 
joan eh (aS) | fds 
oe 
- 
As a result, the element of volume will acquire an angular acceleration 
0 > i 4 49 \= 
oS such that 27é hP edi =falane nla) se Since the motion is a 
2 


| ae —_ a) of VC 
forced oscillation with frequency a the substitution © = (o] S 
ce 


can be made, and the equation of motion becomes 


do 4 . de 7 , 
Ae? Ta dz “eles =O 


e 
The solution of this equation is Oz [ad \ | where J 





; : = 3 
and Y are Bessel functions of the first order, A and B are arbitrary 


Y2 


constants to be determined by the boundary conditions and CL is -g. v2? \ 
vu 
. At the inside wall of the cylinder containing the liquid, 
| Jue 
the motion can be represented by Q@(R> Op since the liquid is 
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following the forced motion of the cylinder which has amplitude GS, - 
the suspended lucite cylinder experiences a shearing torque equal to 
ce ry oe due to the liquid and a restoring torque — Oy due to 
the Scan wire which has a torsional constant K. As a result the 


bob acquires an angular acceleration such that 


all WY? (33},.’ moO, = 71 OC, 


where I is the moment of inertia of the cylinder. Putting the boundary 
conditions into the solution for one can arrive at the equation of 


motion for the inner cylinder. 


a Gar? [F@HY ad - FH) 4G] 
(i tor) LY, Gy) TR) —Tery;, (4 @) abi? (FOR) Yau-Tan Voi) 


where |= “gar hn and TOY is The Yalye of Ay, (a2) acZir 
(tw -h) d (a2) 


Although this represents the exact relationship between the desired 
; es: Se) ee : , 
quantity and the measured quantities 1t 1S in an inconvenient 
OR 
form since { , the unknown, is buried in the complex argument of the 


Bessel functions. Making a Taylor expansion of the functions with 


argument 4 about Q ly the following result is obtained: 


es: 
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and neglecting higher order terms results in 


OK — ln] ett = et? [cos +4 Sing} 
Sener) <a t-4) ye : 
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a COS y + Sing] = | ae 9) | ail : 
in [COSY +g sina] 4 [Gi 48Pj)w - & } ots | (moet) <p | 


42 now may be found by the solution of a complex quadratic equation or, 
if the term involving yy amy be neglected, then by a first order equa- 
tion. 


p= GoSP -j eng ah (Gite -G) =O 
yy rm 7 h ( Take 3) 


In this investigation it was possible to neglect all terms of 


l 
higher order than ~ , because of the errors in the measurement of 


phase angle and amplitude. 
All of the constants may be calculated from direct measurements 
except for the torsional constant of the suspension wire. The torsional 


constant may be determined from the frequency of oscillation of the 


inner cylinder suspended in air using the relation 


ae 
bh 2 (2 Hf \ on 
The values of the constants for the three inner cylinders used 


are shown in Table l. 
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4. Experimental Procedures 

Basically, the operation consisted of filling the inner cylinder 
with the sample of Polyox to be tested, then immersing and centering 
the inner cylinder. The relative magnitude of oscillation and phase 
lag angle is then recorded as a function of oscillation frequency. 

The samples of Polyox were all mixed in a 2000 ML beaker with a 
three bladed beater at 100 RPM. The powdered Polyox was slowly added 
during the mixing cycle and care was taken not to form clumps in the 
solution. 

The turntable was carefully leveled before a run to insure that 
the inner cylinder would be suspended concentrically with the outer 
cylinder. 

The inner cylinder was centered before the Polyox was introduced 
by checking its position over a pointed plug which fitted in the outer 
cylinder center. The inner cylinder then was lowered to a prescribed 
depth after the outer cylinder had been filled. 

The response of the elastometer is very dependent on the viscosity 
of the fluid, the distance of Separation between the cylinders, and the 
frequency of oscillation. As an example, in a 1% solution of Polyox 
with a 12 mm. diameter inner cylinder and a 14.4 mm. diameter outer 
cylinder at frequencies of less than 3 Hertz, both cylinders appear 
to move together, and at frequencies higher than 8 Hertz, the motion 
of the inner cylinder becomes very small. Accurate measurements of 
viscosity are difficult under either of these conditions. The radius 
of the inner cylinder, or the restoring force, must be varied in order 
to cover the frequency range of interest. In this investigation, dif- 


ferent inner cylinders were used with the same outer cylinder. 
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Most of the problems encountered were mechanical in nature, such 
as the centering of the inner cylinder, and the stability of motor 
operation. 

Although the apparatus was constructed with a constant tempera- 
ture jacket, this was not employed during the runs. The measurements 
were made at room temperature, which was 24 eT oCr = ine dynamic vis- 
cosity of Polyox is temperature dependent but Chester's results indi- 
cate that any uncertainty caused by a + 1° temperature variation would 


be less than the experimental error. 


17 


5. Results and Discussion 

The results are best shown graphically (Figures 4,5) with the 

\A 7 Fold 

plots of (t, the dynamic viscosity, and ‘§ the modulus of rigidity, 
vs. frequency for different concentrations of Polyox. Selected values 
are also presented in Table 2. These two quantities are plotted with 
a straight line approximation for different concentrations. The dynamic 
viscosity decreases with increasing frequency of oscillation while the 
modulus of rigidity increases. 

; / 

ge and G are also plotted (Figure 6) as a function of concen- 

; ini (ee 
tration at a frequency of 3 Hertz. Both {V, and change greatly with 
increasing concentration and both of the plots have the same general 
shape. 

An exact estimation of the amount of error is difficult to make in 
this investigation. The elastometer was tested with a Newtonian solu- 
tion of known viscosity and the results checked within 10%. The error 
appeared to be random, so it could not be used as a corrective factor 
to be applied to other measurements. 

: 7 

The most critical variable was the phase lag angle. Since ie 

7 ° e e ° e 
and G are direct functions in the computation formula of the sine 
and cosine of the phase lag angle, they were markedly affected by any 
errors in that angle. This was particularly true when the angle 
measured was near O) I or 7! where the sine or cosine change most 
2 
rapidly. 

The phase angle is determined by measuring the ratio of thé sweep 
delay to the period. For small delays the uncertainty in measurement 
due to lack of trace stability and eye estimation is of the same order 
of magnitude as the observed quantity. An estimate of the uncertainty 


in the phase angle is = 3°. 
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A phase angle meter was also used to try to improve the accuracy 
of the measurement; however, it was not capable of giving good results 
on the irregular pulse outputs from the photocells. 

The amplitude radio was easily measured to within an accuracy of 
* 5% since both light beams swung through an are ranging from 1% to 8 
inches in length. 

In a given series of measurements with a particular inner cylin- 
der, the amplitude ratio starts at unity for low frequencies, increases 
to a maximum at system resonance, then decreases as frequency increases. 
The phase lag angle starts at zero, increases to at resonance, and 
approaches as the frequency increases. A value of vy and e was 
determined at whatever frequency the amplitude ratio‘and phase lag 
angle was measured, and these values were used to plot the dependence 
on frequency. An alternate method would be to plot the phaSe lag 
angle and amplitude ratio as a function of frequency and use values 
determined from this plot to calculate 1. and, - Regardless of the 
method used, each point on the plot would have a different estimate of 
error due to uncertainties in determination of frequency, of phase lag, 
of amplitude ratio, and of the other physical constants of the system. 
The points with the largest estimated error are the end points ina 
particular series of meaSurements, i1.e., those at the low and high 
frequencies. The phase lag angle for these points is nearO and [[f 
and any uncertainty in it affects the results greatly. 

The straight lines were plotted by judgment and an overall error 
of * 10% is considered reasonable, based on the results of measure- 
ments of fluids of known viscosity. 


While the investigation was in progress, it was learned that 


Ie) 


Michael C. Williams at the University of Wisconsin had made similar 
measurements on 1.46% Polyox solution using a Birnboim concentric 
cylindrical pumping viscometer. His measurements covered a wider 
frequency band but his results match the results of this investigation 
and those of Lt. Kinner and Lt. Reister very closely. His results are 
plotted with the results of this investigation. 

The simplest way of accounting for the mechanical properties of 
Polyox solutions is to introduce a single viscosity \l, and a single 
shear modulus G—@. One method of superimposing these properties is 
the Maxwell model. In the Maxwell model, the fluid is represented by 
a Maxwellian element consisting of a spring (representing the shear 
modulus) and a dashpot (representing the viscosity) in series. If Paes 
is the zero shear rate viscosity, (O is the frequency of the applied 
disturbance, then the Maxwell model for viscoelastic solutions gives 
a formula for the dynamic viscosity 

vi, le 
Penh 


where 


— ai le 
rape 


Cw = the shear modulus at infinite shear. 





/ 

Using values of Vy determined in this investigation at 3 Hertz 
and 6 Hertz to solve the equation of the simple Maxwell model Corayey 
and 1 for various concentrations of Polyox, values “pal tel Ates com- 
puted which are the same order of magnitude as the low shear rate 
viscosity coefficients measured at 6 RPM with a Brookfield viscometer. 
The values are compared with the relaxation times in Table 3. The 
agreement een ae is better at the lower concentrations and these 


results indicate that in the frequency range from 1 to 11 Hertz the 
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behavior of Polyox solutions may be approximated by a single Maxwell 
element. The results of Lt. Kinnier and Lt. Reister indicate that as 
frequency increases, this approximation is no longer valid. 

The elastometer did not prove to be a good experimental device 
for investigating Polyox solutions. The properties of Polyox solutions 
change greatly with frequency and this requires many combinations of 
inner and outer cylinders to cover a limited frequency range. This 
investigation, I feel, has filled a gap in the Polyox spectrum. Fur- 


ther work is needed in smaller concentrations and at higher frequencies. 
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